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Abstract 

Hadronic event shapes, that is, event shapes at hadron colliders, could provide a great 
way to test both standard and non-standard theoretical models. However, they are sig- 
nificantly more complicated than event shapes at e + e~ colliders, involving multiple hard 
directions, multiple channels and multiple color structures. In this paper, hadronic event 
shapes are examined with Soft-Collinear Effective Theory (SCET) by expanding around 
the dijet limit. A simple event shape, threshold thrust, is defined. This observable is 
global and has no free parameters, making it ideal for clarifying how resummation of 
hadronic event shapes can be done in SCET. Threshold thrust is calculated at next-to- 
leading fixed order (NLO) in SCET and resummed to next-to-next-to-leading logarithmic 
accuracy (NNLL). The scale-dependent parts of the soft function are shown to agree with 
what is expected from general observations, and the factorization formula is explicitly 
shown to be renormalization group invariant to 1-loop. Although threshold thrust is 
not itself expected to be phenomenologically interesting, it can be modified into a re- 
lated observable which allows the jet px distribution to be calculated and resummed to 
NNLL+NLO accuracy. As in other processes, one expects resummation to be important 
even for moderate jet momenta due to dynamical threshold enhancement. A general 
discussion of threshold enhancement and non-global logs in hadronic event shapes is also 
included. 



1 Introduction 



Every collision at the Large Hadron Collider (LHC) will involve radiation of quarks and 
gluons. In fact, the vast majority of LHC collisions will involve nothing else (except a few 
photons), until the quarks and gluons hadronize and decay. These pure QCD events provide 
a critical way both to test the standard and to search for new physics. At e + e~ colliders, an 
excellent way to characterize QCD events is with event shapes. Event shapes are observables 
which are both infrared safe and simple enough that a trustworthy perturbation expansion 
can be performed. They usually involve a single scale, which guarantees that only one type 
of large logarithm can appear. Resumming these large logarithms can then be used to get 
very accurate predictions. For example, thrust was calculated at NNLO [I] then resummed 
to N 3 LL accuracy [2] (based on [21 H]) and fit for power corrections [5] using Soft-Collinear 
Effective Theory (SCET). Together, these results provided some of the best tests of QCD, 
and important new physics constraints [6j. It is therefore natural to attempt the calculation 
of similar observables at hadron colliders. 

At hadron colliders, pure QCD events are significantly messier than at e + e~ colliders. As 
with e + e~, it is natural to look at event shapes which are dominated by kinematic configura- 
tions with two outgoing jets. Even then, there are at least five new ingredients which must 
be understood: 1) The energy distributions in the incoming hadrons are non-perturbative, so 
one must integrate over different initial states weighted by parton-distribution functions. 2) 
The process involves four directions of large energy flow (the two protons and the two jets) as 
opposed to two in e + e~; thus, non-trivial angles are involved. 3) The 4-parton configurations 
involve multiple channels (e.g. qq — > qq and gg — > gg). 4) There are multiple color structures 
for each channel which mix. And 5) one cannot get by with simple global observables; at a 
hadron collider one must either cut out or suppress the contribution from the region near the 
beam which cannot be measured. Such restrictions can generate non-global logarithms which 
complicate the resummation [8] . 

Dealing with the complication 1), the PDFs, is now well-understood. Following the example 
of Drell-Yan, we can integrate over the PDFs in an effective theory. An important lesson 
learned from the Drell-Yan example is that the matching scales in the effective theory should 
be chosen after the convolution with the PDFs rather than before. This avoids issues involving 
the Landau pole singularity of QCD and allows for dynamical enhancement of the partonic 
threshold contribution to the Drell-Yan mass [9]. Similar enhancements have also been seen 
in more complicated signals, such as direct photon [10] and threshold tt production [TT] . 

Complication 2), multiple collinear directions, has been addressed using traditional thresh- 
old resummation in [12] and for the hard matching coefficients in SCET [131 EH E3 HS1 E3 EH 
HP3 I2D] • In order to understand how angular dependence cancels in a physical cross section in 
SCET, a process one step simpler than dijets was studied in JTU], direct photon production. 
Direct photon involves three collinear directions. The SCET factorization theorem for this 
case was derived in [10] , which required a non-trivial cancellation of various angular factors, ti 
production [11] also has multiple directions which must cancel, but the directions are associ- 
ated with the top quark velocities in the heavy quark limit rather than with massless partons. 
In both cases, the cancellation was demonstrated with explicit 1-loop calculations. Direct 
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photon and tt also involve PDFs and dynamical threshold enhancement as well as elements of 
complication 3), multiple channels. 

Complication 4), multiple color structures which mix, was studied at next-to-leading order 
using traditional threshold resummation in [211 [22]. It was a ls° studied in SCET for off-shell 
Green's functions in [131 113 EE3 EEE] and for threshold tt production in [TT]. The dijet case is 
more complicated than threshold tt because of the jet functions, the myriad channels including 
some with identical particles, and the necessary phase space restrictions for the soft function. 
In SCET language, the color mixing occurs only in the evolution the hard and soft functions, 
with the jet function and PDF evolution being diagonal in color space. The hard function 
evolution is universal, and observable independent. The 1-loop Wilson coefficients for all the 
2 — > 2 partonic channels and their renormalization group evolution equations to 2-loops were 
recently calculated in SCET in [19], based on results of [231 EH EE]. The soft function is a 
cross section for emission from Eikonal Wilson lines. These emission graphs depend critically 
on the observable of interest. Therefore, there is a non-trivial check on the consistency of the 
factorization theorem that color evolution of the soft function cancels that of the hard function 
and its renormalization scale dependence is compensated exactly when the hard, jet and PDF 
evolution equations are combined. A general formulation of these constraints was presented 
in [TU] , although no explicit soft function was given. A main result of this paper is to confirm 
that the general results of [19J actually hold through explicit calculation of a specific case. 

Complication 5), the observable, is a difficult one. The types of observables that will 
be most interesting to compute first are probably analogs of thrust. For event kinematics 
in which the final state looks like a pair of almost massless jets, thrust r can be calculated 
from the sum of the jet masses up to higher order terms r = {m\ + m jJ + " • • > where 
mj 1 and mj 2 are the masses of the two jets. Jet masses are important because the jet mass 
distribution is singular at tree-level da ~ 5(m 2 j)5(mj^), and therefore very sensitive to QCD 
radiation. Getting the theoretical prediction for the shape of jet mass distributions to agree 
with data without Monte Carlo tuning probably requires at least the NLO distributions and 
resummation beyond NLL, as was demonstrated with thrust [21 IS] and heavy jet mass [22] in 
the simpler e + e~ case. 

To calculate jet masses, or a generalization of thrust, at hadron colliders one must deal 
with the complication that most of the energy in a typical event disappears down the beam 
pipe and is unmeasurable. Simply cutting out the beam is unlikely to work, as it will generate 
non-global logarithms which cannot be resummed. There are many ways to proceed. In a 
comprehensive study of hadronic event shapes [21] with the program CAESAR [27J , a number of 
suggestions were put forward. For example, global transverse thrust [28J uses only transverse 
momentum, so that particles in the beam direction do not contribute. Some authors have 
suggested using beam functions [21]. Ref. |26j summarizes some of the relevant issues about 
observables, non-global logs and the experimental measurements. 

In this paper, we will take a somewhat different approach, inspired by the practicality 
of dynamical threshold enhancement for the cases of Drell-Yan and direct photon. In these 
cases, one calculates an observable in the machine threshold limit: the protons collide into the 
final states of interest (a lepton pair for Drell-Yan, or a photon and a jet for direct photon) 
plus only soft radiation. That is, there is no outgoing collinear field in the direction of the 
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beam. Although the physical regime of interest looks nothing like this - there is a beam 
with a lot of energy - resummation still seems to be quantitatively extremely useful in the 
physical regime when extrapolated away from the machine threshold. This phenomenological 
observation is partially understood: the logarithms of the partonic threshold variable which 
should be small well away from the machine threshold are enhanced by a factor related to the 
die-off of the PDFs near x — > 1 j9]. To be clear, we will not attempt to demonstrate threshold 
enhancement for hadronic event shapes here. Instead, we simply observe that the threshold 
region is physically motivated, and use it to clarify some of the relevant issues in SCET. 

With these motivations, we will begin in this paper by studying an observable we call 
threshold thrust, which is simply the hadronic generalization of thrust in e + e~ events. It can 
be defined exactly as regular thrust, but the factorization theorem will only apply near the 
machine threshold, where the protons annihilate into two jets and soft radiation only. The 
threshold thrust distribution we define and calculate in this paper is observable, and should 
agree with experiment near threshold. Unfortunately, there will not be any data anywhere 
close to the machine threshold, due to the die off of the PDFs near x — > 1. Thus, the 
observable will need to be modified so that the large logarithms we resum will still dominate 
in the physical regime. Nevertheless, threshold thrust demonstrates a number of new features 
relevant for any hadronic event shape and lets us check the general expression for the color 
mixing RGE presented in [19]. Since threshold thrust has no parameters, such as a jet size 
R, it can be thought of as an immaculate toy model, perfect for studying factorization in a 
purely hadronic environment. 

After establishing and checking that threshold thrust has the expected properties, we 
propose a related observable, asymmetric thrust, which is more likely to be important phe- 
nomenologically. Near threshold, threshold thrust can be computed, up to power corrections, 
by separating the event into two hemispheres and adding the hemisphere masses. Asymmetric 
thrust is defined by, instead of taking two equally-sized hemisphere jets, we take one jet to have 
size R, and the other jet to be everything else in the event. Threshold thrust is not expected 
to be useful away from threshold because when the beam remnants have large energy, their 
contribution will significantly affect hemisphere masses. For asymmetric thrust, the beam 
remnants both contribute to the same jet, and so their contributions largely cancel. Asym- 
metric thrust also satisfies a number of properties which we expect observables undergoing 
dynamical threshold enhancement to have, as described in Section HI 

The organization of this paper is as follows. In Section [21 we set up the calculation 
of hadronic events shapes in SCET. A number of features are universal and apply to any 
hadronic dijet event shape. These universal features include the 1-loop Wilson coefficients for 
the hard function their renormalization-group evolution equations. We also summarize known 
results on inclusive jet functions and PDFs in SCET language, since they apply to many 
observables. Finally, the parts of the soft function RGE which are off-diagonal in color space 
are universal, since they must cancel the hard evolution, so we include a discussion of these as 
well. Next, we introduce threshold thrust, in Section [3j Section [3731 calculates explicitly the 1- 
loop soft function for this observable. While soft function integrals are observable dependent, 
the color factors for emission from an Eikonal line are universal. The way the color factors 
combine with the momentum-integrals to reproduce the color structures in the hard function 
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is non-trivial. So, this calculation demonstrates the kinds of calculations we expect to occur 
in other hadronic event shapes. In Section I3.5[ we combine the ingredients into a closed-form 
expression for threshold thrust which can be evaluated numerically. In Section H] we discuss 
how to go away from the machine threshold. We review previous examples and catalog some 
lessons learned, extracting some general principles for threshold enhancement. This leads to 
the proposal of asymmetric thrust, which is only briefly mentioned, and will be the subject of 
future work. Conclusions are in Section 



In this section, we describe elements of hadronic event shape calculations which are universal, 
postponing until Section |3] an actual event shape definition and event-shape specific results. 
This section is essentially a summary of results derived in [19] and can be skipped if the reader 
is already familiar with the contents of that paper. 

In Soft-Collinear Effective Theory, a cross section is calculated by expanding around a 
threshold limit. We define hadronic dijet event shapes as observables which force that, at 
the exact threshold, the incoming and outgoing partons have massless 4-momenta and that 
there is no phase space left for soft radiation. At threshold, only massless 2—7-2 scattering 
in QCD contributes, which includes the qq — » qq, qq — > gg and gg — > gg channels and 
their various crossings. Going slightly away from threshold allows the partons to become jets 
and for soft radiation to be produced. When the jets have masses much smaller than their 
energies mj <C Ej and the soft radiation has energy of order -E so ft ~ m'j/Ej, the cross section 
factorizes into a convolution of different components which are calculable separately in the 
effective theory. For these event shapes, the cross section will factorize into a form which 
generically looks like 



Here, Hjj is the hard function, Su is a soft function, and J and / are, respectively, jet 
functions and parton distribution functions (PDFs). dll is the differential phase space and 
iVinit reminds us to average over initial states. / and J are color indices which only affect the 
hard and soft functions, since the jet functions and PDFs are color diagonal. All the functions 
in this equation have an implicit channel index, which we suppress for clarity. In the threshold 
limit, one can show using SCET that the channels do not interfere and the cross section for 
each channel can simply be added together [ID] . 

2.1 Wilson coefficients and the hard function 

The first step in an SCET calculation is to match QCD to SCET. This means calculating 
Wilson coefficients for SCET operators so that amplitudes in QCD are reproduced at the scale 
fj, order-by-order in perturbation theory. For example, one of the SCET operators relevant for 
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gg -» qq is 

off- = (x^ilAAlxi) (^D^VV^n)' • (2) 

Here, Pi = |(1 + 75) projects out the left-handed quarks while the gluon helicities are chosen 
explicitly to be + and — and so this operator describes the helicity subprocess g + g~ — > 
QlQl- The part of the operator in the left brackets comprises quark jets \n and gluon jets 
A n ±, which are collinear quarks and gluons plus associated collinear Wilson lines. The right 
part of this operator comprises soft Wilson lines Y n and y n , in the fundamental and adjoint 
representations, respectively. The soft part is independent of spin. i,j are color indices for 
the fundamental representation of SU (3) and a, b are adjoint color indices. For this gg — > 
qq channel, there are 3 color structures and 8 spin choices, although many of the Wilson 
coefficients for these 24 operators are the same. There are additional operators for the crossed 
channels, such as qg — » qg, whose Wilson coefficients can be computed from those of gg — > qq 
via crossing relations. For qq' — > qq' and its crossings there are 2 color structures, and for 
99 ~^ 99 there are 8 color structures. Special care is required in computing the Wilson 
coefficients for qq' — > qq' when q and q' are identical particles. All the channels and crossings 
are cataloged and their Wilson coefficients computed to 1-loop in Ref . [19] . 

The RG evolution of the Wilson coefficients is known to 2-loop order. It has the relatively 
simple form 



d -c T M 



7cus P ^- In — + 1h ) Ojj + 7 cusp M/j 



c T M- (3) 



dln/i 

Here, 7 cusp = 24 + • • • is the cusp anomalous dimension, ch is the hard group Casimir 

c H = n q C F + n g C A , (4) 

where n q is the number of quarks or antiquarks and n g is the number of gluons, including both 
initial and final partons (e.g. ch = 2Cf + 2Ca for gg — > qq). The hard anomalous dimension 
Jh can be calculated as a series in ot s . It is expressible in terms of quark and gluon anomalous 
dimensions as 

1h = n q -i q + n g -i g , (5) 

with 7^ = (j 2 -) (— 3CV) + • • • and 7 g = f^) (—A)) + • • • known up to order [IB] . The matrix 
Mjj describes the color mixing. It depends on the channel and the scattering kinematics 
through logarithms of ratios of the parton momenta, however, it does not depend on a s . 
General expressions for Mjj are given in the color space formalism in [181 US]- The color 
factors and crossings are worked out explicitly in [19]. For example, for the qq' — > qq' channel, 
the mixing matrix is 

'«7,(to=s+i7r)-CU(lii*f + 2i,r) 2(\n=s+in) 
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Explicit results for all 2 — > 2 processes in QCD are presented in [T9"]!*! 

In SCET, the Wilson coefficients from matching to QCD can be combined into a hard 
function by summing over spins. The hard function in any particular channel is defined by 



H 



I j 



r* 
j > 



(7) 



where T labels the spins and / and J label the color structures. An important consequence of 
the color mixing being universal is the existence of a natural basis in which the evolution of 
the Wilson coefficients, and therefore the hard function as well, are diagonal. In this basis, 



din// 



7cus P c H In 



+ X K + \* K > + 27 



YH 



2/9 (a a 



where Xk are the eigenvalues of M/j, and we have used the fact that 7# is real. This equation 
is solved by 



H K K'{s,t,U,fi) 



a s (//) 2 



exp 



2c H S(/i h ,n) - 2A H (fj, h , fj) 



x exp 



A r (//ft,//) ( X K (s,t,u) + X* K ,( s ,t,u) + c H ln 



t 


)] 







H K K'{s,t,u,fx h ) , (9) 



where 



J a s (v) 



/3(a) 



Q d« / 



Ar(v, yu) 



(10) 



and Ah{v,ij) is the same as A r (v, //) but with 7^ replacing 7 cusp . Closed-form expressions for 
these functions in renormalizat ion-group improved perturbation theory can be found in 



2.2 Jet functions 

Next, let us turn to the jet functions. Jet functions come from matrix elements of collinear 
fields in SCET. The simplest jet function is the inclusive jet function, J(m 2 , //) which depends 
only on the mass of the jet. This well-travelled jet function has had phenomenological ap- 
plication to B decays [5Tj . deep- inelastic scattering [52], event shapes [2], and direct photon 
production [10]. It is all we will need for the applications in this paper. The inclusive quark 
jet function can be written to all orders as [33] 

^P 2 ,^) = exp[-4:C F S(fi j ,fi) + 2Aj q (fi j ,fi)}j q (d r]jq )—[—\ — — -. (11) 

P / 1 {Vjq ) 

1 The mixing matrices given in [19) differ from previous results of Kidonakis et. al [21) only by convention- 
dependent diagonal terms and a slightly different basis in the gg — > qq channels. In addition, while the results 
of [21| are derived by computing the virtual contributions to the renormalization of soft Wilson lines with a 
collinear subtraction, the Mjj matrices in [19] are calculated simply from virtual graphs in full QCD. The 
results of |19) also include the finite terms in the matching calculation and the generalization of the mixing to 
NNLO, which is necessary for NNLL resummation. 
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The function j g (L), which is the Laplace transform of the jet function, can be expanded at the 
scale fij order-by-order in a s . Each order is a finite polynomial in L. This function is known 
exactly up to 2-loops, and its complete L dependence to 3- loops. For example, to order a s , 



J 2 / cy 2 

4C F 3C F L + C F 7 

2 I 3 



(12) 



The placeholder rjj in Eq. ffTTl) is to be evaluated at r)j = 2CpA r (fij, /x) after the derivatives 
are taken, with A^(y, fi) given in Eq. (fTUj) . Aj (i/,fi) is an evolution kernel like Ap(z/, /x) and 
Ah(v,h) above. The gluon jet function has a similar form, with different coefficients [31]. 
Details of all of these functions and an example application can be found in [TO] . 

To check the renormalization group equations for threshold thrust below, it is easier to use 
the RGE satisfied by the jet functions, rather than to use its solution. The RGE, for either 
quark or gluon jets, is 



dln/i 



ji (Q 2 ,v) 



Q 2 

-2Cjj.7 CUS p In — - 27 



(13) 



where C Rq = C F , C Rg 
in [TO]. To order a,, 



Ca and the anomalous dimensions and j% functions can be found 



47T 



-3C F ) + ■ 



(14) 
(15) 



2.3 Parton distribution functions 

In any theoretical calculation at a hadron collider, parton distribution functions /i/jv(^, A*) for 
parton i in nucleon iV will play some role. These PDFs are non-perturbative, but their renor- 
malization group evolution equations, the DGLAP equations, are perturbative and known to 
3-loops. In general, the evolution mixes all the different parton species. However, near the 
endpoint x — > 1, species mixing is suppressed and the evolution equations simplify. Observ- 
ables such as threshold thrust which are calculated in the threshold region where x ~ 1 allow 
us to use this simplification to prove renormalization-group invariance. In practice, the full 
non-perturbative PDFs with the general x evolution can be used for phenomenology, with the 
x < 1 evolution compensated for at fixed order by careful matching. Again, see [10] for more 
details and some quantitative results. 

We define the Laplace transform of the parton distribution functions by 



fi/N (C A 1 ) 



dx exp 



1 — x 



fi/ N (x,ji) . 



The evolution equations near x = 1 are then 

d 



dln/i 



frniC, y) = [2C Ki 7 cusp in c + 27^] f i/N ((, fj) . 



(16) 



(17) 
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Expressions for jf q = jf q and 7 Js to 3-loop order can be found in (TO]- To order a s the PDF 
anomalous dimension are simply negative of the a s jet-function anomalous dimensions 

^ = (3C F ) + ■■■ (18) 



Air, 

7/9 = (£) m + ■ ■ ■ • (19) 

However, at higher orders, the PDF and jet function anomalous dimensions are independent. 
2.4 Soft function: general observations 

The soft function in a hadronic event shape calculation depends on the definition of the 
event shape and must be calculated for each observable separately. However, a number of 
features of these soft functions are universal, and will apply for any observable. Because the 
renormalization group evolution is color diagonal for the jet functions and PDFs, the color- 
mixing terms in the soft function evolution must exactly compensate the color-mixing terms in 
the hard function evolution. This was explained in more detail in [TJ5]. Here, we review those 
general results for completeness. A direct calculation of the threshold thrust soft function is 
given in the next section, which may help clarify the notation introduced here. 

The soft function is a matrix in color space. It is calculated from matrix elements of Wj 
which are time-ordered products of the soft Wilson lines appearing the the SCET operators: 

Sij ({*}, <) = Yl (°l \ X °) ( X °\ |0> F s ({k}) , (20) 

where the sum is over soft radiation in the final state. The function Fg{{k}) encodes the 
dependence on various projections on the soft momenta related to the definition of the ob- 
servable. An explicit example for F$ will be given when we define the observable. No matter 
what the observable is, the soft function can only depend on directions of the various Wil- 
son lines, and on arbitrary soft scales relevant to the projections. Because of factorization, it 
cannot depend on the energy of the jets, the hard scales s, t, u or the energy fractions Xj of 
the PDFs. 

As was shown in [19], based on insights in [181 EI]> f° r the factorization theorem to hold, 
the soft function must satisfy 



d 



dln/i 



Su({Q},n?,n) = -^(TO^^^rL-rfX^TO^f^) , (21) 



where 



rfj = ^7cus P cq In + 7cu Sp ?"«) + Is^j $u + 7cuspM 7 j(nf ) . (22) 

This soft function RGE has much in common with the RGE for the Wilson coefficients, Eq. (J3]). 
In particular, Mjj is the same and in both cases all the color mixing is proportional to 7 CU sp- 
The Casimir cq controlling the Sudakov logs and the remainder function r(nf) depend on 
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the channel but are independent of a s . The soft anomalous dimension 75 depends on the 
observable and the channel but not on the color structure. For any particular observable, the 
factorization theorem will fix 75 to be a linear combination of hard, jet, and PDF anomalous 
dimensions. 

The general solution to the soft function RGE in Laplace space is 



- 2c Q S(fi s , n) + 2A s (fi s , fi) 

S K K>({Q}X,t*»), ( 23 ) 

where Xk are the same eigenvalues of Mjj used for the hard function evolution. This can be 
transformed back to momentum space using the techniques described, for example, in [2| [18]. 

We will now proceed to apply the general results of this section to a particular example: we 
define a hadronic event shape, work out its factorization formula, calculate the soft function 
at order a s , check renormalization group invariance, and produce a closed form expression 
resummed to NNLL. 



Skk>{{Q},K,h) = exp 



x exp 



A r (fis, H) + X* K , + r«) + r«)* + 2c Q In 



W± 
11. 



3 Threshold Thrust 

In the previous sections, we presented results applicable for resummation of any observable 
which is dominated by dijet configurations. We saw that the color mixing terms in the renor- 
malization group evolution can be diagonalized in a closed form to all orders in perturbation 
theory. Now we will proceed to work out the details for a simple observable, threshold thrust. 
Some advantages of this observable are that it involves only inclusive jet functions and that 
all radiation contributes, making it manifestly free of non-global logs. 

The factorization theorem for threshold thrust is a straightforward combination of the 
ingredients used for direct photon factorization |10] and the e + e~ thrust factorization theorem 
in SCET [3J. For simplicity, we will present only the physical derivation rather than the 
technical one. 



3.1 Kinematics 

Thrust in e + e~ is defined by 

r = l-max^M^, (24) 

where the maximum is taken over directions n. This somewhat complicated definition obscures 
the fact that when an event has two final state jets, r is simply the sum of the masses of the 
jets normalized to the machine energy, up to power corrections of higher order in those masses. 

In the threshold limit for a hadronic collision, we have two incoming protons of momenta 
Pi and P2 annihilating into two jets with momentum P£ and P£ and soft radiation. In the 
threshold limit, there are no beam remnants and the jets are back-to-back. Since the are no 
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beam remnants, the final state looks just like a possible final state from e + e~, and so thrust 
can, in principle, be measured experimentally without modification. For simplicity, we define 
the jet momenta as the vector sum of the momenta of all the radiation in the left and right 
hemispheres, respectively. The hemispheres are defined with respect to the thrust axis, the 
maximal n in Eq. (124"]) . which becomes equal to the jet directions in the threshold limit. Which 
hemisphere is called L or R is arbitrary. To further simplify, rather than using the definition 
in Eq. (|24|) , we will simply define our threshold thrust variable r as the sum of the hemisphere 
masses, normalized to the machine energy, r = [Pi + Pj|] . 

Our labelling convention is that we will use Pf and P% for the incoming protons and p^ 
and P2 for the incoming partons. We call the outgoing hemisphere momenta, and P£, 
while the parton level quarks and gluons will be P3 and p%. We allow for either p% or p^ to 
align with either Pj* or P£. The hadronic Mandelstam variables are 

S=(P 1 + P 2 ) 2 , T = (P 1 -P R ) 2 , U = (P 1 -P L ) 2 , (25) 

and the partonic level Mandelstam variables are 

s = (pi + p 2 ) 2 , t = 0i - p 3 ) 2 } u = (pi- Pi) 2 , (26) 

with the usual definition of the momentum fractions 

rf = x 1 P?, p% = x 2 P£, (27) 

At leading order, the partons are all massless. We can then define lightlike 4-vectors in the 
direction of the 4-momenta, so 

p^ ~ E x n%, p% ~ E 2 n%, p'i, ~ E 3 n$, p% ~ P 4 < , (28) 

where E 1 ,E 2 ,E i , and E± are the energies of partons. 

At threshold E x = E 2 = E 3 = E 4 = VS/2 = \fs/2 = E C m, Til = n 2 and n 3 = rT 4 . In 
this limit, there is only one dimensionless Lorentz-invariant ratio on which the hard and soft 
functions can depend: 

ni-ns n 2 -riA , . 

v = — — - = -^—^ > . 29 
ni-n.4 n 2 -n 3 

In terms of s, t, u, 

u = -, l + u = —. (30) 
u —u 

These relations will be necessary to check the factorization theorem. 

To understand threshold thrust, it is helpful to pursue the analogy with direct photon 
production [10] . For direct photon, the threshold observable was M\ = S + T + U, where M x 
is the mass of everything-but-the-photon. In that case, Mx can be written as a function of 
only the photon p T and rapidity y, M\ = E^ M — 2p T E CM coshy. For dijets, we can also look 
at S + T + U but now we find that 

S 4 = S + T + U = P| + Pi. (31) 
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So, in both cases S + T + U gives the observable of interest. In the photon case, it is the 
photon pt and rapidity. In the dijet case, it the sum of the hemisphere masses, which is equal 
to threshold thrust times the machine energy squared: 

r=§ = ^ + T + t/] . (32) 

In fact, if we write this in terms of the transverse momentum and the rapidity y of the jets, 
we find an expression identical to direct photon: 

r = l ^=coshy. (33) 

So, by calculating threshold thrust, we will produce the px and rapidity spectrum of the 
jet. Note that the two jets always have the same pt and opposite rapidities even away from 
threshold since they are hemisphere jets. 

Since the kinematics are like those of direct photon production we will express the final 
distribution in terms of the variables v,w,pr and y, which are related to s,t,u,x\ and x 2 
by pEJEED], 

s = -, t = , u = — — , p T = — (34) 

WW WV V s 

Pt 1 y Pt 1 _„ , QK v 

II = 7f^ e ' l2 = vl* e ' (35) 

with v = 1 — v. We will use these definitions implicitly in the following. 

To write down the threshold thrust distribution in SCET, we now use the fact that thresh- 
old thrust is just e + e~ thrust integrated over appropriate parton luminosities with colored 
initial states. So we define the partonic thrust variable 

S4 = m R + rrij = s + t + u = — , (3o) 

V w 

This is equivalent to the partonic variable m 2 x used for direct photon in [10]. Then we just 
write down the e + e~ thrust distribution, convoluted with PDFs0 

7T~ l ^ dV l dW ~2 [ X lfl/Ni(.Xl, fi)] [X2f2/N a (x 2 , fi)] 

iV init JSTeV J^T.eV W 



dp T dy 8np T . ,„„ 

channels ^fs y/s v 



X 



dm\dm 2 4 dk J{m\, //) J{m 2 A) [i)SKK'{k, v, /i)5(s 4 — m 2 , — ml — a/s/c) (37) 



2 The ^2 factor in the integrand is a convention for power corrections. It is chosen so that if the integrand 
were written as ds 4 instead of dw, all the nonsingular S4 dependence is in the PDFs. This is slightly different 
from the convention chosen in [10 j which had all the nonsignular w-dependence in the PDFs. 
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In contrast to e + e~ thrust, the hard and soft functions now have color indices, because the 
initial states are colored. They also can depend on s,t and u. In the threshold limit, where 
the hard and soft functions are calculated, there is only one independent dimensionless ratio 
of these quantities, which we have taken to be v. 



3.2 Threshold limit 

The threshold thrust distribution calculated with SCET is only formally valid as r — > 0, 
equivalently, as s — > S. It is in this limit that the hard and soft functions are calculated and 
the factorization formula can be checked. Studying this limit will also clarify the definition of 
the soft function and the way the different hard directions interact. 

To derive the factorization formula, we need to relate the physical quantities, defined at the 
hadron level, to perturbative quantities, defined at the parton level. The momenta (1 — xi)Pf 
and (1 — x^jP^ represent the momenta from the two protons which are not involved in the hard 
interaction. In the threshold limit, this initial-state radiation is soft. We can then define left 
and right so that that the remnants of the first proton, (1 — #i)Pf go into the left hemisphere 
and the remnants of the second proton, (1 — X2)P 2 go into the right hemisphere!! We also know 
that in the threshold limit, the hemisphere momentum scales like a collinear field and looks 
like a jet, i.e. its mass is much smaller than its energy v/Pf <C Er- So, each hemisphere must 
contain initial-state radiation from one of the proton remnants and final-state radiation from 
one of the jets. This is shown graphically in Figure [TJ We see that there are two cases, when 
u > t, the left hemisphere aligns with the p% direction and when t > u the left hemisphere 
aligns with the direction. It follows that, in the threshold limit, we can decompose the 
total hemisphere momentum as 

P£ = (P^ + k$)G(t -u) + {P? + K)Q(u - t) + (1 - x 2 )P$ (38) 
p£ = (pf + fc£)e(i -u) + (if + fc£)e(u - t) + (1 - Xl )P? . (39) 

The sum P^ + denotes the allocation of final-state radiation into collinear and soft sectors. 
This separation is not well-defined and k% and k% must be integrated over to form a physical 
observable. 

The right hemisphere mass is then 



p2 



ml + 2E 3 (n 3 -k 3 )-u(l-x 2 ) + --- 6(t - u) (40) 



+ 



mj + 2E 4 (n 4 -h)-t(l-x 2 ) + --- Q(u-t). (41) 



and similarity for P|. The ellipses denote terms of higher order in the ki or (1 — Xj). We have 
written to, = a/ (if M ) 2 to remind us that these are jet masses. Since only one component 
of the soft momentum in each hemisphere contributes, we we will abbreviate the projections 



3 We could instead have defined left and right so left is aligned with parton 3 and right with parton 4. 
However, associating left and right with the proton directions is more obviously physical - it does not require 
the experiment to distinguish which parton is which. The threshold thrust distribution is the same either way. 
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u > t 



t > u 



Figure 1: The hemispheres are defined so that the remnants of proton 1 go left and the 
remnants of proton 2 to right. For u > t, the 3-jet is in the right hemisphere and the 4-jet in 
the left hemisphere, for t < u the jets switch sides. 



with k 3 = (n 3 -k 3 ) and /c 4 = (n 4 -/c 4 ). Also, in the partonic center-of-mass frame, the energies 
of the jets are 2E 3 = 2E4 = ^/s. Putting all the ingredients together, we get a cross section 
of the general form of Eq. (PQ) : 



da 



dPMP? 



T -i- 



i,J,r, 
channels 



dxida^dmndm^dfcadfci 



x 5 
x5 



x (C L j* ■ Sj^h, h) -C\ )J 3 (mi)J 4 (mi)f 1 (x 1 )f 2 (x 2 ) 
ml + \fs{n 3 ■ ks) — u(l — x 2 ) Q(t — u) + m\ + A/s(n 4 • /c 4 ) — t(l — x 2 ) ®(u — t) — P\ 



m 



\ + \fs{n 4 ■ fc 4 ) — u(l — X\) <d(t — u) + m\ + \/s(n 3 ■ k 3 ) — t(l — Xi) <d(u — t) — P 1 



(42) 



Here, the sum is over color structures / and J and spins Y. Cj are the Wilson coefficients for 
matching to SCET, determined by virtual graphs in full QCD. Ji{m 2 ) are quark or gluon jet 
functions, depending on the process, and fi(xi) are quark or gluon PDFs. These objects were 
all introduced in Section [2j 

Finally, threshold thrust is the sum of the hemisphere masses. So, 



d 2 CT 



R L 1 dP^dPl 



5 T 



r R + L 
S 
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— — / dxida^dmgdm^d/c (43) 



channels 

x Hj j S n ( k) J 3 (mj ) J 4 (m^ ) f ± (xi ) / 2 (x 2 ) 

x 5 (7713 + ml + y/sk — min(t, u)(l — X\) — min(t, u)(l — x 2 ) — St) , 
where the threshold thrust soft function is defined as 

S^ T {k) = J dhdhSj^h, k 4 )S(k - k 3 - h) . (44) 

We will now show how to calculate all of the objects in this formula and check the renormal- 
ization scale independence to order a s . 



3.3 Soft function 

Now that we have defined an observable, threshold thrust, which is the sum of the hemisphere 
masses in the threshold limit, we can define the soft function more precisely. In this case, it 

Si j(h, h) = (0| W] \X S ) (X s \ Wj |0> 5 (n 3 • - k 3 ) 6 (n 4 • P?° - h) , (45) 

x 3 

where P^ s and P* s are the sum of the momenta in the particles in state X s which go into 
the hemispheres containing parton 3 or 4 respectively. 

As an example, consider the qq' — > qq' channel. For this channel the W/ are constructed 
out of soft Wilson lines Y n in the fundamental representation. The two color structures are 

Wi = T{(y 4 W 2 )(y3V°yi)} 

w 2 = T{(y 1 t iy 2 )(y 3 t iy 1 )}. (46) 

Here, T{} stands for time ordering and 1 is the identity operator in SU(3). The W/ have 
external fundamental indices, which have been suppressed for conciseness. These indices get 
contracted with each other when the WV are combined into the gauge invariant soft function. 
The more explicit expression is 



X 



S u = <0|T {(Y-MY 3 ) l h.AY}T\Y,) %} 

'(Y}TjY 2 ) i4 i2 (Y^TjYi) * 3 n } |o) 6 (n 3 -P*° - k 3 ) 6 (n 4 -Pf s - h) , (47) 



x (X 



where T x ® T x = r a ® r a and T 2 ® T 2 = 1 (8) 1. 

The hemisphere soft function is complicated by the necessary projections into the hemi- 
spheres, which prevents us from simply summing over the intermediate states. At tree level, 
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these projections are trivial, the Wilson lines are numbers, Y n = 1, and the soft function for 
qq' — > qq' evaluates to 



- 1 t::!S t ;: 1 T :'::!it; 1 1 



Tr[r 6 ]Tr[r 6 ] Tr[l]Tr[l] 



2 U 

o c\ 



The color factors here are related to the normalization of the Wilson line in Eq. (1461) . Careful 
monitoring of the normalization is required for the factorization theorem to check. 



3.3.1 NLO soft function 

At order a s , there are contributions to Su from virtual and real emission graphs, however, in 
dimensional regularization, the virtual graphs for matrix elements involving only soft Wilson 
lines are scaleless and vanish. The real emission contribution involves emissions from any of 
the Wilson lines and absorption into any other. The necessary diagrams can be drawn as cuts, 
as shown in Figure |5J The calculation can be split into two parts: calculation of the integrals, 
which we call Is, and calculation of the color factors, which we call Du. The 0(a s ) result 
can then be written as 

Sf^°(k 3 ,k 4 ) = J2 I s(n a ,n b ,k 3} h)D IJ (a } b) . (49) 

a, b 

The sum is over assignments of a and b to any of the four Wilson lines 1,2,3,4. The color 
factors are matrices in color space and depend on the channel being considered. The integrals 
depend only on kinematic factors not depend on the color nor channel, i.e. whether the lines 
are quarks, anti-quarks or gluon. A similar computation was done for ti production in [TTj . 

The integral for emission from leg leg and absorption into the leg, in d = 4 — 2e 
dimensions, is 

Is(n a ,n b M,h) = 9l (^) / (^ e ^)^^) K n ; )( ^. g) 

x |e(n 3 -g - n 4 -q)6(k 4 - n 4 -q)5(k 3 ) + 6(n 4 -g - n 3 -q)6(k 3 - n 3 -q)5(h)\ . (50) 

The second line denotes the projections into the two hemispheres: its first term says "If the 
component of q 11 going to the 4 direction is larger than the component going to the 3 direction, 
then the soft radiation must be going to the 4 direction. It therefore only contributes to k^, 
by an amount equal to n 4 • q." The second term is for the case where the radiation goes in the 
3 direction. At order a s the radiation can either go in the 3 or 4 direction, but not both. 

Because of the explicit appearance of rig and n^, the integral has a different form for 
the various assignments of n a and n b to the various Wilson line directions. To simplify the 
expressions, we can use the fact that rig and 77,4 are back-to-back as are and n^. It is also 
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1 




Figure 2: Diagrams contributing to the soft function at NLO. 

true that the integrals are Lorentz invariant as well as invariant under separate rescalings of 
and n^. Using these observations, the integrals can only be functions of the single ratio 
v — n\ -n-zjnx -ri4, as in Eq. f l2"§j) . The integrals are not smooth functions of v at v = 1, due 
to the convention chosen for the color basis; a different convention might move the poles from 
the t- to u- channel in a particular color operator. We give all the integrals for all values of 
v in Appendix |A] Some of the cases already exist in the literature. For example, if n a = n$ 
and rib — n 4, then the integral is the same one from the hemisphere soft function from thrust, 
which was calculated in [4] and [36]. Another case can be checked against results in [37]. In 
all cases, our results agree with previous work. 

In addition to performing the integrals, the color factors Djj(a,b) have to be worked 
out. These color factors come from traces of products of group generators coming from the 
expansion of the various Wilson lines to NLO. They depend on which line is emitting and 
which is absorbing, and therefore on the assignments of a and b. The NLO color factors for 
all channels are given in Appendix |A] The color factors depend only on which Wilson lines 
are doing the emitting, not on the details of the observable (unlike the integrals). Therefore, 
these color factors should be generally useful for the calculation of any dijet hadronic event 
shape, not just threshold thrust. 

To compute the soft function for a crossed process one must permute which Wilson lines go 
into the color factors with respect to Wilson lines used for emissions in the associated integral. 
To do this, we write 

Djj(o,6)-)-£) JJ (x(a) > x(6)). (51) 

In general, there are 4! = 24 permutations, but they are not all independent. For the qq — > qq 
there are 6 independent permutations xi a ) which are listed in Table [TJ The other 18 channels 
are equal to one of these using either charge conjugation (q <-> q for all 4 quarks) or q •<->■ q 1 . 
The permutations for the gg — > qq channels are listed in Table [2J In this case, there 12 
different permutations corresponding to the 12 different physical channels. For some channels, 
like qq — > gg, there is an ambiguity as to whether it is an stu channel with x( a ) = 4321 or an 



16 



12 ->■ 34 



cross 



X{a) 



12 ->■ 34 



cross 



qq' — > qq' 



sttt 



1234 



2C F \n 



min(t,u) 



gg' —7- g'g 



sut 



1243 



2C F ln 



min(/j,u) 



gg gg 
(gg gg) 



uts 



1432 



2CVln 



min(/j,«) 



gg' — > q'q 



tus 



1423 



min(t,u) 



qq — )■ g'g' 
(gg -> gg) 



1324 



2C F (ln 



- min(t,u) 



gg g'g' 



ust 



1342 



2C F (ln- 



- m'm(t,u) 



— ?7T 



-27T 



Table 1: Crossing relations for the color factors in the soft function for the gg — > qq chan- 
nels. Identical particle channels are bracketed. The "cross" column show the permutation 
of Mandelstam invariants used to cross the Wilson coefficients. The xi a ) columns are the 
permutations used in Du{x{ a )i x{b))- The r(nf) columns show the explicit crossed form of 
r(nf) appearing in Tfj, in Eq. (1581) . 



sut channel with x( a ) = 4312. We can pick either convention, and as long as the appropriate 
stu or sut hard function is used, the physical cross section will be the same. For gg — > gg 
there is only one channel, so it does not get a table. 

To be clear about the x{ a ) notation, consider two examples. For the gg' — > qq' channel, 
corresponding to stu and x( a ) = 1234 in Table HJ the soft function is 



S 



1234 
IJ 



/s(m, n 2 )D u (l, 2) + I s ( ni , n 3 )D u (l, 3) + I s {n u n 4 )D u (l, 4) 



+ I s (n 2 , n 3 )D u (2, 3) + I s {n 2 , n 4 ) J D/j(2, 4) + I s {n 3 , n 4 )D u (3, 4) 



(52) 



where the 2 comes from Du(a,b) = Djj(b,a). For gg — > q'q', which according to Table [T] is 
tsu with x{ a ) = 1324, the soft function is 



QtSU 

°7J 



I s (n 1 ,n 2 )D IJ (l, 3) + I s (n u n 3 )D u (l, 2) + I s (n u 7i 4 )Djj(l, 4) 



+ I s (n 2 , n 3 )D u (S, 2) + I s (n 2 , n 4 )D/j(3, 4) + I s (n 3 , n 4 )D/j(2, 4) 



(53) 



and so on for the other channels. In our notation, q' and q are different flavors of quarks 
(for example, up and down). The soft function does not care if the particles are identical. 
However, for identical particles, one must choose a convention for which color basis to use. 
Our convention agrees with the convention in [19] , and the appropriate permutations are listed 
in Table [U 

Since the soft functions are calculated with real emissions integrals, they are all real. Thus, 
the crossings cannot be computed by expressing the final soft function in terms of stu and 
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12 34 
99 -»■ 99 

qq^gg 



cross 



1234 
4321 



r(nf) 



c m In 



t 



2 111 -i 



12 34 

-> <?<? 

QQ -> 



cross 



1243 
3421 



c m In ■ 



-> £9 
#9 -> 92 



3214 
1432 



C, In 



min(i,«) 



-»■ qg 
gq -»■ #9 



tits 



3241 
1423 



Cm In 



min(t,it) 



-> 92 

gq^ gq 



tsu 



4231 
1324 



c m (In 



min(t,u) 
— Z7r) 



qg -> 
qq^qg 



ust 



4213 
1342 



c m (In 



— min(t,u) 
271") 



Table 2: Crossing relations for the color factors in the soft function for the gg — > qq channels, 
with columns as in Table ED c m and cq are given in Eqs. (160|) and ( 16T1) . and cq has been 
dropped when it vanishes. There are more soft functions for gg — > qq than for qq — > qq since 
Du(l, 2) ^ L>/j(3, 4) in this case. 



then crossing stu, as we do for the Wilson coefficients. In fact, there is no obvious relationship 
between the soft functions for the different crossings and they simply have to be computed 
separately by permuting the arguments of Du(a, b). 

As an example, combining the integrals and color factors for the qq' — > qq' channel, the 
threshold thrust soft function is 



Su(k) = S(k) 



\C f Ca 



a 




4CVln 



(i+^) 2 



8C 2 F C A \nuQ(l 



-8C \C f ln(l + v) 



-%C A C F \n{l + v) 
WC\C F hxvQ{v- 1) 



(54) 



Here, [l/A;]^'^ is a star distribution (see jl]). For the /i-independent part, cfj(z/), see Eq. f ll34jl . 
Other channels will have a different Cjj(u). Also, for other channels, there may be an additional 
piece proportional to [In k/k]^'^. It just so happens that the coefficient of this additional term 
vanishes for qq' — > gg' (and its crossings). 

There is already a non-trivial check we can make from just this expression. The angular 
variable v is related to the Mandelstam invariants by Eq. ( 1301) . At intermediate states of 
the calculation terms of the form ln(u — 1) appear (see for example Eq. (11021) ). These terms 
cannot be expressed simply as ratios of Mandelstam invariants, as is required to cancel the 
scale dependence of the hard function. The check is that all the ln(v — 1) terms drop out 
in the final expression, which involves intricate cancellations among different directions and 
different color factors. 
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3.3.2 Soft RGE 



The RGE for the threshold thrust soft function must be of the general form in Eq. (12 lj) . To 
see that this equation is satisfied, we first transform to Laplace space via 

Su(Q,fi) = J dkexp (-^A Sjj(k^) . (55) 

For example, for the qq' — > qq' channel, this gives 



/ or, \ ( AC F In g - 8C 2 F C A In i6(« - t) 8C A C F In ^ \ Q 

Utt/ I 8C,CVln^ 16C 2 ,C P ln^eft- uW 11 /i ' 1 J 



8C A C F m^ 16ClC F ln^e(t-u) / A* 



We have used Eq. (|30|) to write the soft function in terms of s,t and w. We can now check 
that the soft function for qq' — > qq' satisfies the RGE, Eq. (I2T!) with 

r?j = p.\ ( - c^( ln I + - 8C > ln S (« - 8 ( ln V + ™) \ (57) 



4tt/ V ^F(\ n ^L + i n ) 8C F ln^9(t- 



This matrix is not uniquely defined by the NLO soft function, In fact, the soft RGE would also 
be satisfied for a Tjj with somewhat different off-diagonal terms or imaginary parts. However, 
it is a powerful check on the soft function that it can be written in this way, which is of the 
form expected by RG invariance, Eq. (I2"2"j) . 

For all channels, we find that the soft functions for threshold thrust satisfy the general 
RGE, Eq. fl2~Tj) . which takes the particular form 

rfj = ^7cus P cq ln (^j + 7cuspKnO + 7s ^ 8u + 7cus P M/j , 



with 



r{nf) 




ln ^ — min(t, u) j 


Cm 


= E C * 






initial 




C Q 








initial 


final 


Is 


= + 0(a 2 s ) 


) 



CQ 

2 ; 



(58) 

(59) 
(60) 

(61) 
(62) 

with L(x) = ln(— x) — i7fQ(x), so L(t) = ln(— t), as in [19]. For the crossings of Tfj, direct 
calculation shows that only the t in the L(t) part of r(nf) and the s,t and u in M/j get 
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crossed. The min(t, u) and s factors are the same for all channels. To avoid any confusion, 
r(nf) for all crossings is given explicitly in Tables [T] and [2J The Mjj matrices are the same as 
for the hard evolution, which are given for all channels in [T5] , 

The appearance of min(t, u) in Tfj foreshadows the satisfaction of the RGE consistency 
check we perform in the next section - it will cancel the factor of min(t, u) in the factorization 
formula, Eq. fj43|) . The soft function min(t, u) factor comes from differences in the locations of 
phase space singularities when t > u and when t < u, which trace back to our convention for 
how the fermions are contracted in color space the original operators. There can be no factor 
of min(i, u) in the hard function, since the hard function is analytic and independent of the 
observable. 

To solve the RGE, we diagonalize Mjj. Consistent with [19], we use the index K for 
the diagonal basis. The combinations of Wilson lines for which the soft function evolution is 
diagonal are then Wk, and we write 



W K = F KI Wi 



(63) 



and call the eigenvalues A a-. Since the same off diagonal terms appear in the hard and soft 
evolution, the same basis simultaneously diagonalizes the evolution equation of the hard and 
soft functions: 



Hrk 1 = (F • H ■ F^) KK i S KK i = [(F x ]kk' 



(64) 



The eigenvalues, A# and matrices Fm which diagonalize Mjj are given for all channels in |19j . 
For example, for qq' — > qq', 



W± = X ± Wi + 



C A 



In h iTT 

s 



and the eigenvalues are given by 



W 2 



(65) 



C A _ u 1 











\ -t 




In h %tx 


* 


In h itt 




In h itt 


+ 


s 




s 




s 





C iln 2 - 



(66) 



In matrix notation, we have W± = F±j ■ W/ where F is 

F±i = 



A + g[ln^ + m] 



(67) 



The operators and coefficients for the other channels and color structures are presented in 
Appendix [A] 

Once the soft function is diagonal, its RGE can be solved in Laplace space, as in Eq. f f23|) . 
In this case, since the soft function has a single scale, we write 



S KKI ({Q} } n^ } n) = S KK/ (\n—,u,ij,) 



(68) 
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Then, 

S KK ,(k, v,\i) = exp - 2cqS{h s ,h) + 2A s (fi s ,n) + A r (fJ. a ,fi) (X K + \k> 

t 



(69) 



x exp 



A T (ns,v) 2c m ln 



min(t, u) 



cq In 



where rj s = 2CqA r (/j ls , //). For the crossed processes, only t gets crossed, while the s and 
min(t, u) factor are s and min(t, u) for all channels. Other examples of this notation can be 
found in [21 1 



3.4 Checking RGE invariance at NLO 

Now that we have computed all the ingredients, we can combine them together to form the 
final distribution. Before accounting for all the kinematic factors, we can check that the 
factorization formula, Eq. ( H3l) . is RG invariant at order a s . This must hold for each channel 
and each spin separately. 

To check RG invariance, we first go to Laplace space. We define the Laplace transform of 
the cross section as 

da [°° , f tS \ da 
Then RG invariance requires, for each channel and spin, that in the threshold limit, 



d ln/i 



5> = 



,K,K' 



Q 2 



min(t, u) 



H/2 



Q 2 



min(t, u) 



Q 2 



x H U (^)S K 'K 



Plugging in each of the RG equations, this implies 



0. (71) 



Tcusp 



2ci In ^r—r, — 7 + 2c 2 In — - — 2c 3 In ^— — 2c 4 In ^ 



min(t, u) 



min(i, u) 



+ c H In -* + X K 4- Ai, - 2 (c, In + c m In — L-j + 3 In -L) - A, - A;, 

+ 2 7fl + 2 7/2 - 2 7Jj - 2 7J4 + 2 7H - 2 7s = . (72) 

To check this equation, we need to use that 

ch = c\ + c 2 + c 3 + c 4 (73) 
cq = Ci + c 2 - c 3 - c 4 (74) 
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C m — C\+ C2 



(75) 



These relations show that the part of Eq. (1721) proportional to 7 CU sp is satisfied. The remainder 
is satisfied if 

75 = lh + 7/ 2 - 7J ;i - 7j 4 + 1h ■ (76) 

Using the anomalous dimensions we have calculated, this holds at order a s . This equation 
can the be used to determine the soft anomalous dimension to a?. 



3.5 Threshold Thrust distribution 

The final threshold thrust distribution is achieved by combining the hard, jet and soft functions 
and adding appropriate kinematic factors, as in Eq. (1371) . The result is 



- £ 



1 



da 

dp T dy %irp T ^ N Mt Jp^ eV 

channels y/s 



dv 



!>%ley 

Vs f 



V 

dw— [xifyn^X!, (],)] [x 2 f 2 /N 2 (x2, V>)] 
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Here, C3 and C4 are the fundamental Casimirs (Cp or Ca) for the 3 or 4 jet and 
V = V3 + + Vs = 2(c 3 + c^Arifij, /1) + 2c Q A r (fi s , //) 

We have also used 

I min(t, u) I = p T ^~, — rr 
mm(f , v) 

I v — - = - to write things directly in terms of t> . For the different channels, one must use 
the appropriately crossed hard, jet and soft functions and PDFs, and additionally cross the 
explicit factors of t. 
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4 Dynamical Threshold Enhancement 



The threshold thrust variable we have calculated in this paper allows us to calculate the 
Pt and rapidity distribution of dijet events jets near the machine threshold. In this threshold 
region, the final state consists of two almost massless jets and nothing else. Such a threshold is 
physically irrelevant because the final state in any realistic collision has beam remnants which 
contribute to the hemisphere masses, preventing them from ever being close to massless. 
However, since the beam remnants have no transverse momentum, the pt distribution of 
realistic jets should have singular contributions very similar to the hemisphere jets in threshold 
thrust. Thus, it should be possible to calculate the physical jet pt spectrum using SCET and 
to compare directly to data. 

The reason the SCET factorization theorem is applicable only in the threshold region is 
because in that region the momenta of the incoming partons approach the momenta of the 
incoming hadrons. In this limit a hadron can be thought of as a parton with the other, 
spectator, partons in the hadron providing only power suppressed contributions. Thus, the 
singular logarithmic terms are guaranteed to be the dominant contribution to the cross section 
at threshold. For the large logs, and their resummation, to be important away from threshold, 
their relative contribution must be somehow enhanced from what one might naively expect. 
This is called dynamical threshold enhancement. 

Dynamical threshold enhancement has been shown to hold for other processes such as 
Drell-Yan and direct photon production. In Drell-Yan, the observable is the mass of the 
final state leptons, = {pf+ + p^) 2 and the threshold region is when — > S. In Direct 
photon, the observable is the mass of everything-but-the-photon, M\ = S + T + U and the 
threshold region is when M\ — > 0. Dynamical threshold enhancement in direct photon can be 
understood intuitively. Mx includes contributions from the mass of the final state jet. Away 
from threshold, the factorization theorem does not guarantee that the jet mass would be small 
compared to its energy. However, we know from observation that the typical jet mass in a 
realistic event is in fact small. Therefore, the large logarithms at threshold related to the jet 
mass are also relevant in the physical regime. 

For threshold thrust, the observable r is very similar to the observable for direct photon, 
since r = (S + T + U)/S. We expect large logarithms of the jet mass to be similarly important 
even away from threshold, since typical jet masses are small compared to their energy in 
both cases. However, there is an important difference. For direct photon, we can write the 
final state momentum as pg nal = p^ + p x . Since p x is the momentum of everything but the 
photon, it includes both beam remnants. Thus, the direction of p x is the same at the parton 
and hadron level, even though its energy may be vastly different when s < S. That the 
direction is the same is critical for SCET, because the direction is a label for the operator, 
and appears in the final distribution. In contrast, the energy of the final state jet cancels out 
at an intermediate stage of the calculation, and only reappears through the matching scales. 

For dijet events, we can write the final state as pg nal = Pr+Pl- Since one beam remnant 
goes into one hemisphere and the other beam remnant goes into the other hemisphere, the 
direction of the final states at the parton and hadron levels can be vastly different. In other 
words, the hemisphere axis is sensitive to the boost of the event, while the direction of the 
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jet in direct photon is not. If the direction label on the jet changes, the matching coefficients 
should change too. So, there is no reason to expect the threshold thrust calculation to be 
relevant in the physical regime. 

Suppose we try to modify the definition of the observable so that its direction is stable 
away from threshold. For example, we could incorporate a jet definition, say a cone of size 
R. Then instead of threshold thrust, which is the sum of hemisphere masses-squared, we 
look at the sum of the squares of the masses of the two hardest jets of size R in the event, 
T/2 = {Ph) 2 + (Pl) 2 - Since these jets would not include the beam, they should be stable between 
the parton and hadron levels. This modified observable has (at least) two problems. First, 
it is not global, so there may be relevant non-global logs which we are not resumming, even 
at threshold. Second, the vanishing of the observable does not guarantee the dijet kinematics 
relevant to the factorization theorem. While the observable vanishing does imply that both 
jet masses must vanish, it does not enforce that the remaining radiation should be soft. For 
example, there could be other hard jets in the event. Thus, while the observable is good away 
from threshold, it is not useful at threshold. 

In summary, there appear to be two general principles which an observable should satisfy 
for dynamical threshold enhancement to work. We would like to have a single definition of 
the observable which can be applied both at threshold and away from threshold so that 

1. The direction of the final state jets is the same at the parton and hadron levels. 

2. The vanishing of the observable implies the kinematics associated with the factorization 
theorem. 

In the case of hadronic dijet event shapes, the second point means that only 2 — » 2 processes 
should contribute when the observable is exactly at threshold. For Drell-Yan, both of these 
criteria hold: there are no final state jets, and when m\ — > S, the remaining radiation must be 
soft. For direct photon, these criteria hold: the jet defined as the everything-but-the-photon 
is always back-to-back with the photon, and Mx — > forces the jet to be massless and the 
remaining radiation to be soft. Threshold thrust satisfies the second condition but fails the 
first because the beam remnants change the jet direction. 

In addition, to avoid complications involving non-global logs 

3. The observable should be global. 

This means all hadrons in the event should contribute to the observable [7J [8] . For direct pho- 
ton and Drell-Yan, the observables are the photon 4-momentum or lepton pair 4-momentum 
respectively, which seem clearly not to be global since they do not involve the hadrons at all. 
However, in direct photon, what is actually calculated in SCET is the mass of everything-but- 
the-photon, Mx, which is global. In Drell-Yan, one can also rephrase the calculation as being 
of the mass of everything but the lepton pair, M\ = S + T + U — m 2 u which is global as well. 
For threshold thrust, M\ = S + T + U = Pj, + Pf , which is also global. This is not to say 
that observables which are non-global must have non-global logarithms, or that non-global 
logarithms cannot be resummed. Moreover, there are additional subtleties within the class of 
global observables, such as whether they are directly or continuously global [8J. So it makes 
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sense to start in SCET by studying event shapes which are as global as possible, so that all 
the large logarithms will be threshold logs. 

Now let us consider some modifications of threshold thrust. First, take the example 
mentioned above, the sum of the squares of the masses of the two hardest jets of size R, 
Tj2 = (pi) 2 + (Pr) 2 - It will satisfy condition [1] but not condition [2j Even though the jets are 
forced to be massless when rj2 = 0, the radiation outside the jets is not forced to vanish by 
Tj2 = alone. The variable r + = j;(p R + p1) 2 , discussed in [TJ] and [3H] has the the opposite 
problem. In this case, r + = 1 does force the radiation outside the jets to vanish, however, 
the jets are not forced to be massless. They can have masses as large as rrij/Ej < R. Al- 
though this observable may be interesting, resumming logs of the jet masses cannot be done 
without also resumming logs of R. If R is as large as a hemisphere, then r+ = 1 exactly, so 
there is nothing to calculate. A related observable, suggested by p2] and studied in [39] is 
r_ = j;(2p R -p L ) = r + — Tj2. When this observable is at threshold, t_ = 1, the final state 
is forced to have two massless jets and no out of jet radiation. It satisfies both criteria [1] 
and [2J and therefore we expect it to undergo dynamical threshold enhancement. However, it 
is not global: because the jets have size R, only the radiation within the jets contributes to 
the value of Pr-Pl- Unlike direct photon or Drell-Yan, the mass of everything not in the jets 
is also non-global beacuse it excludes the hadrons in the jets. Thus, when attempting NNLL 
resummation, there may be non-global logs in dcr/dr_ of the same order as the logs which are 
resummed. 

Finally, consider the following modification of threshold thrust. Instead of taking two 
hemispheres, one finds the hardest jet of size R (using whatever algorithm one wants). Then 
the radiation in the event is divided into radiation in the jet, which is summed into a mo- 
mentum 4- vector p? n and out of the jet, which is summed into p% nt . Then one computes the 
asymmetric thrust 

TA = \ Mnf + (rfut) 2 ] • (80) 

This will allow us to resum logarithms of the pr of the jet, which can be compared to data 
and should undergo threshold enhancement. The exact .R-dependence can be included through 
careful matching to a fixed order calculation. 

For R the size of a hemisphere, this observable is just threshold thrust. For R small enough 
so that neither beam is included in the jet, this observable will satisfy condition [TJ In fact, 
if R is calculated as a distance in rj — <p space, the beams will never be included in a jet 
since they are at rj = ±oo. Condition [2] is satisfied as well. When Mx — > 0, both the jet 
mass \J (Pi n ) 2 and the mass-of-everything-else ■sj (Po Ut ) 2 are forced to vanish. The initial state 
radiation must also be soft, for the same reason as for threshold thrust. Thus, at threshold, 
only 2 — » 2 scattering contributes, as desired. Finally, the observable is global, so there are 
no dangerous non-global logs to worry about. 

5 Conclusions 

In this paper, we have presented the first computation of a hadronic event shape in pure QCD 
events to next-to-next-to-leading logarithmic accuracy. This calculation of this observable, 
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which we call threshold thrust, involves understanding a number of issues which were not 
present in previous work. The hadronic event shapes we have been concerned with here are 
dominated by dijet configurations. In the singular limit, at threshold, only 2—^2 scattering 
contributes, but there are many channels and many color structures in each channel. The 
evolution in color space is complicated, but has a number of simplifying features. In particular, 
there is a color basis in which the evolution is diagonal. The evolution of the hard functions 
was studied previously in [H], and constraints on the RG evolution for the soft function were 
derived. Here, we have computed the soft function for threshold thrust explicitly to 1-loop. 
This both confirms the general expectations from JT9] , and provides the NLO finite parts of 
the soft function necessary for NNLL resummation. 

Although threshold thrust is infrared safe and observable in principle, the large logarithmic 
contributions to it which the SCET calculation provides are unlikely to have a significant effect 
in the physically relevant kinematic region. The threshold where the calculation is valid has 
x 1 for both hadrons, so that the hadron momenta and the hard scattering parton momenta 
approach each other. It is not unusual for the threshold logarithms to be important away from 
threshold. In fact, for related processes, such as Drell-Yan and direct photon production, this 
is known to happen due to dynamical threshold enhancement. We isolate two features of 
observables which are likely to be important for dynamical threshold enhancement to occur: 
1) The jets in the observable should have the same directions at the hadron and parton levels 
and 2) The observable should force the kinematics of the factorization theorem. In addition, 
to avoid non-global logs, 3) The observable should be sensitive to radiation everywhere. We 
defined an alternative to threshold thrust, asymmetric thrust, which satisfies these criteria. 
While threshold thrust is the sum the masses-squared of hemisphere jets, asymmetric thrust 
sums the masses- squared of a jet and everything-but-the-jet. The critical feature of this 
modification is that both beams end up in one side, so their contributions to the threshold 
observable largely cancel. It will be interesting to study the phenomenology of asymmetric 
thrust, which, after matching to the exact NLO distribution, should allow us to compute the 
jet pt spectrum at NNLL+NLO accuracy. 

Many of the results and lessons from this paper will help us understand more general 
hadronic event shapes. While the pt of a jet in dijet events is an important observable, there 
are many other observables which one could compute using the same 2 — » 2 hard kinematics 
which would provide complimentary tests of QCD. Some of the global observables discussed 
in or 2-jettiness which has been proposed using SCET [10], may be calculable at NNLL 
using the results of this paper and of [IS] without too much more work. There are also many 
interesting non-global observables which can be studied and measured using the same dijet 
sample. For example, it would be nice to calculate the mass of a jet in dijet events. Studies 
in e + e~ collisions of thrust [21 [5], which is the sum of jet masses, and the related heavy jet 
mass [25], have shown clearly that resummation at NNLL and beyond is critical to getting 
some shapes correct. Because of the beam remnants, a jet mass distribution at hadron colliders 
requires a less inclusive jet definition than at e + e~ colliders. Much progress has already 
been made studying exclusive jet definitions in SCET [UJ [37J [38], and phenomenological 
applications at hadron colliders are surely just around the corner. However, to go to NNLL 
requires an understanding of not only non-global logs in SCET, but also of super-leading 
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logs [42], which arise as violations of color coherence in processes with at least four hard 
colored particles, such as the dijet configurations we have been studying. It will be interesting 
to see if SCET can provide a consistent framework for resummation and precision calculations 
of general non-global observables at hadron colliders. 
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A Threshold thrust soft function 

The calculation of the soft function involves taking matrix elements of Wilson lines. To order 
a s in dimensional regularization only the real emission diagrams contribute. These can be 
drawn as cut graphs, as in Figure [2j There are contributions where the gluon comes out of 
the n a leg and is absorbed into the n b leg, where n a and n b are any of ni,ri2,n 3 or 714. Since 
the graphs with n a = n b vanish, there are 6 possibilities. The calculation can be split into 
two parts: calculation of the integrals, which we call Is, and calculation of the color factors, 
which we call Djj. The result can then be written as 

Su(k R , k L ) = Df?5(k 3 )6(h) + J *K, n b , h, k 4 )Djj( X (a), X (6)) + 0(a 2 s ) . (81) 

a,b 

The notation x( a ) an d x(^) indicates to permute the a and b indices in the color factor, Du, 
with respect to the integral, Is, for the appropriate crossing. The crossings and permutations 
are listed in Tables [I] and [21 We work in d = 4—2e dimensions with MS subtraction throughout. 

A.l Integrals 

The master integral for the soft function is 

Is(n a ,n b ,k 3 ,h)= 9s ^— j J J^y-^il )% ) (Wa . g)(w6 . g) 

x |e(n 3 -<7 - n i -q)5(k i - n 4 -q)5(k 3 ) + Q(n 4 -q- n 3 -q)5(k 3 - n 3 ■ q)S(h)j , (82) 

where n a and n b are the directions of the two Wilson lines where the gluon attaches, which 
can be any of m, ri2, n 3 or 77,4, as in Figure [2 Thus, there are six different integrals to be 
done. As discussed in the text, the integrals can only depend on one Lorentz invariant ratio 

ni-n 3 n 2 -n 4 

v = = , 83 

ni • n 4 n 2 ■ n 3 

which greatly simplifies the calculations. 
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At order a s the soft gluon can either go into the hemisphere containing Wilson line 3, or 
Wilson line 4. Thus, at order a s the master integral can be written as 



Is{n a ,n b ,k 3 ,h) = 



/i 2£ /i 2e 

S(k 3 )L ab (u) + —^8{k A )R ah (v) 



l+2e 



34) 
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J7) 



These integrals do not depend on £, by dimensional analysis, but some expressions tend to 
appear simpler before rescaling I away. We now proceed to calculate the 6 dimensionless 
integrals L ah {y). 

The integrals are easiest to do in a basis where the 5 and G functions are easy to integrate, 
which means choosing lightlike coordinates along the 3 and 4 axis. So we write 
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In this basis the measure, including the S and G functions, becomes 
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?l-2e 



d d gG(g )(5(g 2 )G(n3-g - n 4 -g)<5(fc 4 - n A -q){- • •} 
where fi d = 27r d / 2 /r(f ) and 



d-3 



dx 



^ r , , s 



sin 2e 



n 3 -g 



and # is the angle between gj_ and n 4 j_. Then, 



-^afe(^) 
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n a -n b 



ix 4Ji fj„ sin 2£ 6 1 (n a • g) (n& • g) 

The inner products that appear in the denominator are, in this basis, 



n\ • q — ix + v — 2\Tvx cos o\ 

H 1 + v I i 



(91) 



(92) 



(93) 
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n 2 -q 



Y^{ X1J + 1 + 2 cos 6>j (94) 

n 3 ■ q = g_ = £x (95) 
n^ - q = q + = £ . (96) 

Since n x = n 2 , the angle between gj_ and n 2 ± is n — 9 thus changing the sign of the cos 9 
term in the expression for n 2 ■ q. The terms 1 + v arise in the center of mass frame where 
n 3 + n 4 = (2, 0) since 

(ni-n 3 ) 2 2 

1 + v = 1 + = ? , = . 97 

(ni-n 4 ) (ni-n 4 ) (n 2 -n 3 ) 

There are 6 ways to choose two of the four directions for n a and rib, with n a ^ rib. These 
integrals are tricky, but can be done in a straightforward manner. Doing the 9 integral first 
leads to a form involving a hypergeometric function. The x integral then may have singularities 
at x = 0, x = 1, x = v and x = oo. Also, when v = 1, the and directions are on the 
hemisphere boundary. Thus, special care is required around v — 1 and it turns out that the 
regions with v > 1 and ^ < 1 have different functional dependence on v. All of the singularities 
are regulated by e, but special care must be taken to properly treat the various branch cuts. 
It is helpful also to use some tricks similar to those described in 



case n a = n,3, rib = n± 

This case is the same as for the e + e~ hemisphere mass distribution [U 136] . 



(e lE 7rY f°° dx 



7T J 1 X l+£ 



Expanding to order e gives 



, s 2 7T 2 

= --e—. 99 

£ O 



case: n a = n l5 = n2 
For this case, 



^12^) = fid-s(l + ^ 

7T 



,2 



f°° /" r d6> x~ £ 

x dx — ~ = = -. (100) 

J i J Q sin 9 (x + f — 2^/ ux cos 9){xu + 1 + 2^/uxcos9) 

This integral is poorly behaved at v = 1, where and are on the hemisphere boundary. 
However, the integral is symmetric in v — > -, so we can write the result as 

L 12 {v) = F 12 (v)e(v - 1) + F 12 (I) 6(1 - !/) , (101) 
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and just perform the integral assuming v > 1. The result is 
F 12 (u) = ~ + 4 In (1 + v) + 2e ^-3Li 2 H-j - 2Li 2 (1 - v) 

571-2 \ 

-2 In 2 (z/ + 1) - In (u - 1) In u + 2 In (u + 1) In v + 2 \n 2 (u) - — ) ■ (102) 
Note that Li^iv) is real for all v. 
case: n a = m, 712, = 113,714 

This is the most complicated case and we will start with n a = n\. The integrals are most 
easily evaluated using the same light cone basis as in previous integral. Then, 

fe 7B 7r)" [°°, v r d9 

3 / ax- / -r-^ n , . — 

71 Ji x Jo sin 9 [x + v — 2 v vx cos 1 



LvM = ^^-3 / dx U - / ; - - n m (103) 



M") = L T 1 ^/ da; / „:, 2e ^„ I .. o /r^^m - (1° 4 ) 



(e 7B 7r) e o /- 00 r d0 

I2 d _ 3 / dx / —^7^ , — 

^ A Jo sm (x + v — 2 v z/x cos 1 

The results are written in terms of four functions in order to account for the separate cases 
v > 1 or v < 1. We find 

LxM = F 13 (y) Q(y - 1) + G 13 (y) 6(1 - u) (105) 
L u (v) = F u (u) Q(u - 1) + G u (u) 6(1 - v) , (106) 



where 



F 13 (v) = ~ + 2 In (v) + 2 In (1/ - 1) 

+2e ^-3Li 2 Q - ln 2 (z/ - 1) + 2 ln(u - 1) ln(i/) - 2 ln 2 (z/) + ^ 

G 13 {v) = -2 In (1 - v) + 2e ^ ln 2 (l -u)+2 ln(l - 1/) \n(u 

FuM = 21n + 2, (-Li 2 (1) - In 2 " 1 

2 „ , , „ ( 2 ,. ^\ 



G 1A {y) = - - 2 In (1 - v) + 2e \Li 3 (u) + ln 2 (l - U )-—J • (107) 
The case n a = n 2 can be expressed in terms of the same integral as the case n a = n\. 

L 2i {v) = L u {^) = G u e{v - 1) + F u (^J 6(1 - !/) (108) 

£ 23 H = L 13 (i) = C?i 3 A) 6(1/ - 1) + F 13 6(1 - v) . (109) 
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A. 2 color factors 



The color factors that accompany each integral depend on the channel and the crossing. They 
are particular to the basis of operators used and the associated color structures in the soft 
functions. We only need to calculate the color factors Djj(a,b) for the three representa- 
tive channels (qq' — > qq',gg — > qq and gg — > gg), with the crossed processes determined by 
permutations of the indices, as described in the text and in Tables [T] and [2J 

case qq' — > qq' 

The general definition of the soft function was given in Eq. ( I4"5j) . For the qq' — > qq' channel, 
the Wilson lines are 




(110) 



with Ti = r a and T 2 = 1. A more explicit form of these operators, Eq. fj47|) is 



S u = J2(0\t{(YMYs) h iMT}Y,) \} \X S ) 



x (X S |T {{YlTjY 2 ) S^TjY,) \} |o)5(n 4 -P 4 x -fc 4 )5(n3-P3 X -fc 3 ) . (Ill) 



The tree-level color factor is 




(112) 



and the 1-loop color factors are 
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case gg — > qq 



For the gg — > qq channel, the Wilson line structure is 



w? = {ylT^y'Yi = F 3 f (yfT/ b 'yty y 4 



(116) 



where T" 



i are 
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T ab = T b T a 
T ab = S ab ^ 



(117) 



and is a soft Wilson line in the adjoint representation. The explicit form for the soft 
function for this channel is 



yl{ylT\y^ y 3 

(ylTjy 2 y b y 4 



ba 



x ( X 



'2 



(118) 



) 5 (n 4 -P 4 x - h) 5 (n 3 -P 3 x - k 3 ) 



where we have written the color indices explicitly, with a, b being adjoint indices and ii,i 2 
being fundamental indices. The tree level color factor is 
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The 1-loop color factors are 
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D JJ (1,4) = D JJ (2,3) 
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case: gg gg 

The Wilson line structure for gg — > gg is 
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where the T abcd are given by 



rpabcd 


= Tr [r a T b T c r d } 


rpabcd 
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= Ti[T a T c T b r d ] 
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1 2 


= Tr [rVW] 
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The soft function for this channel is 
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The tree level color factor for = 3 is 
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The 1-loop results, also with N = 3, are 

/ 



D JJ (1,2) = D JJ (3,4) 



12 



20 


2 


65 


-7 


2 


-7 


27 


4 


2 


20 


2 


-7 


65 


-7 





4 


65 


2 


20 


-7 


2 


-7 


27 


4 



\ 

27 


-7 -7 -7 -16 -7 -16 -27 -6 -27 
2 65 2 -7 20 -7 48 27 
-7 -7 -7 -16 -7 -16 -27 -6 -27 
27 27 -27 -27 18 -18 
48 48 48 -6 48 -6 18 144 18 
27 -27 27 -27 -18 18 J 



33 



D JJ (1,3) = D JJ (2,4) = ^ 



D I j(l,4)=D IJ (2,3) = ^ 
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A. 3 Final threshold thrust soft function 

With the integrals done and the color factors calculated, all that remains is to combine them 
together to form the threshold thrust soft function. The threshold thrust soft function is 
determined from the hemisphere soft function by 



Su{k) = / d^d^S/jOs, k 4 )8(k - k 3 - fc 4 ) 



a,b 



Rab(v) + L ah (u 



Du(x(a),x(b)) 



(130) 
(131) 



This is straightforward to compute from the results in this Appendix, so we will just give a 
few representative examples. To compute the regulated soft function, we expand in e using 
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(132) 



and drop the - and \ terms using MS subtraction. The is a star distribution (see for 

example [1] for a definition and some relations). 

In the qq' — > qq' channel (stu), for which x(a) = a, the threshold thrust soft function is 
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S u (k) = 5(k) ' 2 



\C f Ca 



o ci 
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4C F ln 
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The constant (/i-independent) part, Cjj(u) are messy functions, which take different form for 
v > 1 and v < 1. For example, the 11 component of the qq' — >■ gg' threshold thrust soft 
function is 



c f 



(C 2 + 2) Li 2 (^j + (C 2 A - 2) ( - InV - 1) + In ^-1) + 4Li 2 (l - z/) 

+ (1 - 3C 2 A ) In 2 v + (AC\ - 6) lnz/ ln(v - 1) + 4 ln 2 (z/ + 1) 
- 41nz> In(i/ + 1) + 7T 2 e(i/ - 1) 

- (2C 2 +1) (Li 2 (l-z/)+lni/ln(l-z/)) + (C 2 A - 4) In 2 i/ + 4 ln> + 1) 



- 41m> + 1) lni/ - ln(l - i/) + y (2C 2 + 7) 

ILi, I ^— ^ +ln 2 (l ) 



0(1-1/). (134) 



The other constants have similarly complicated expressions and we see no reason to write 
them out explicitly 

To compute the soft function for a crossing, one sums the integrals against the color 
factors with the Wilson lines permuted as in Tables [1] and [2j For example, the soft function 
for qq — > q'q' (ust) comes out as 



\C f Ca 



s (/,) ^il^^j + g)^) 



-4C F ln[(l+u)u]+8C F C A In uB{u-X) 
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These Cjj(u) are different from the Cjj(v) for qq' — > qq', but not worth writing out explicitly. 
In the gg — > gg channel, 

= c A cf c A c F S(k) + (^) J 4j{v)5{k) 

\ CaCf CaCf 2C\Cf J y 



~t1 C F 



CaCf 



(c A c% 

I -°f C A C% CaCf ] | 16(C A - C F ) 



In 



k l IM 



16C A lnz/0(l - v) 







k 
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/ 2C F C A \nv+C A \n{l+v) -lav AC A \n{l+v)+2C A In v 

+AC F C A ( C A ln(l+^)-ln<y -2C A In v+4C A ln(l+v) 

\ 2CAlni/+4CAln(l+j/) -2CU In v+iC A Ml+f) 8C^ln(l+y) 

The functions cjj(v) are again different from the other channels, and to messy to write out. 
Note the additional [\n(k)/k}^'^ factor in this channel, which happened to be absent for the 
four quark channels. 
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